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Abstract 


Measurements  of  the  nonlinear  viscoelastic  behavior  of  two  adhe¬ 
sives,  FM-73  and  FM-300,  are  presented  and  discussed.  Analytical  meth¬ 
ods  to  quantify  the  measurements  are  given  and  fitted  into  a  framework 
of  an  accelerated  testing  and  analysis  procedure.  The  single  Integral 
model  used  Is  shown  to  function  well  and  Is  analogous  to  a  time-tempera¬ 
ture  stress-superposition  procedure  (TTSSP).  Advantages  and  disadvan¬ 
tages  of  the  method  accessed  with  the  latter  residing  primarily  with  the 
nature  of  the  creep  power  law  which  is  utilized. 
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Chapter  1 
INTRODUCTION 


Background  Information 

The  connection  or  the  joining  of  components  is  an  integral  part 
of  engineering  structural  design.  This  is  especially  true  of  the 
aerospace  industry,  where  many  thousands  of  structural  components  must 
be  joined  to  form  a  single  aircraft.  Any  decrease  in  the  overall 
structural  weight  of  an  aircraft  results  in  a  corresponding  increase 
in  its  allowable  payload,  thereby  making  the  aircraft  more  cost 
efficient.  This  economic  incentive  has  been  a  major  factor  behind  the 
increased  usage  of  lightweight  metallic  and  composite  materials  in 
structural  aircraft  components  in  recent  years.  For  the  (high  strength/ 
low  weight)  "advanced"  composite  materials,  conventional  fasteners  such 
as  rivets  are  inefficient,  as  they  sever  the  strong  reinforcing  fibers. 
In  the  case  where  metallic  components  are  to  be  fastened,  riveting 
causes  high  stress  concentrations  around  the  rivet  holes,  reducing  the 
fatigue  resistance  of  the  particular  components.  An  additional  draw¬ 
back  of  riveting  is  that  rivets  are  sometimes  made  of  steel  and  there¬ 
fore  may  add  significantly  to  the  total  weight  of  the  structure.  One 
way  to  avoid  the  use  of  rivets  and  their  attendant  problems  is  the 
use  of  adhesive  bonding  to  fasten  structural  components. 

The  desire  to  minimize  weight  as  well  as  the  economy  and  ef¬ 
ficiency  of  construction  in  a  variety  of  industries  has  created  the 
recent  interest  in  high  strength  structural  adhesives.  Structural 


components  are  often  subjected  to  extreme  variations  in  temperature  and 
humidity,  as  well  as  stress.  Obviously,  any  relatively  new  procedure 
such  as  adhesive  bonding  must  first  go  through  an  exhaustive  testing 
and  analysis  program  before  actual  implementation  is  allowed.  One  of 
the  largest  in  depth  research  efforts  conducted  to  better  understand 
adhesive  bonding  was  the  "Primary  Adhesively  Bonded  Structures  Tech¬ 
nology"  (PABST)  program  [1-6]  conducted  by  McDonnell  Douglas  for  the 
USAF  Flight  Dynamics  Laboratory.  The  PABST  program  spanned  a  period 
of  five  years  (1975-1980)  and  has  been  described  as,  "one  of  the  most 
successful  technology  development  programs  ever  undertaken"  [3]. 

Some  of  the  advantages  of  adhesive  bonding  over  riveting  (sum¬ 
marized  from  references  [1-23])  are  as  follows: 

*  Force  distribution  over  a  comparatively  large  area  resulting 
in  relatively  low  stresses 

*  Reduction  in  total  weight  of  the  structure 

*  Possibility  of  eliminating  high  stress  concentrations  at  rivet 
holes 

*  Favorable  vibration  damping 

*  Joint  resistance  to  corrosion  (with  proper  surface  preparation) 

*  Reduced  initial  and  maintenance  costs 

Although  the  future  of  adhesive  bonding  looks  quite  promising, 
there  are  a  few  problems  associated  with  the  analysis  and  design  of 
bonded  joints  which  must  be  solved  before  adhesive  bonding  can  be 
utilized  to  its  fullest  capabilities.  The  most  important  of  these 
problems  is  the  continued  use  of  the  single-lap  bonded  test  coupon  as 


the  basic  standard  for  determining  the  thin  film  shear  properties  of 
adhesively  bonded  joints.  The  single-lap  test  forms  the  basis  of  most 
military  adhesive  specifications  as  well  as  being  a  current  ASTM 
standard  (ASTM  D  1002).  The  results  obtained  via  this  standard  test 
have  been  shown  to  vary  with  the  thickness  of  the  adhesive  layer  as 
well  as  with  the  type  of  adherend  material  used.  An  additional 
peculiarity  of  the  single-lap  test  is  that,  "even  though  it  is 
described  as  a  shear  test,  the  failure  load  is  rarely  if  ever  influenced 
by  the  shear  strength  of  good  structural  adhesives"  [4].  Clearly,  a 
better  standard  test  specimen  is  needed  if  thin  film  adhesive  shear 
properties  are  to  be  determined  accurately. 

In  recent  years  there  has  been  an  increasing  trend  toward  using 
the  thick-adherend  short-overlap  test  coupon  [20]  as  the  basis  for 
determining  accurate  thin  film  shear  stress-strain  adhesive  response. 
With  the  development  of  the  KGR-1  extensometer  [20]  by  R.  B.  Krieger  at 
American  Cyanamid,  accurate  shear  strain  measurement  for  the  thick- 
adherend  test  specimen  is  now  possible.  Although  the  thick-adherend 
short-overlap  test  is  generally  considered  to  be  superior  to  the 
single-lap  test  for  determining  thin  film  adhesive  shear  properties 
[24],  it  is  not  an  official  ASTM  standard  at  the  present  time. 

An  additional  specimen  geometry  presently  under  investigation  is 
the  cracked-lap-shear  (CLS)  specimen  [18,19].  The  CLS  specimen  was 
developed  by  Lockheed  to  simulate  Mode  I/Mode  II  ratios  typical  of 
bonded  aircraft  structures  [19],  The  CLS  test  specimen  is  currently 
being  subjected  to  round  robin  examination  by  several  ASTM  committees. 

At  this  point  it  appears  as  though  the  CLS  specimen  may  turn  out  to  be 


a  very  useful  test  geometry  for  determining  accurate  in  situ  adhesive 
properties. 

In  general,  most  adhesives  are  polymers  and  recently  epoxy  based 
materials  have  become  very  popular  [25-28].  The  nature  of  polymers  is 
such  that  they  have  viscoelastic  or  time  dependent  moduli  and  strength 
properties  which  are  greatly  influenced  by  environmental  conditions, 
especially  temperature  and  humidity.  The  time  dependent  nature  of  the 
adhesive  causes  the  bonded  joint  to  be  time  dependent  as  well.  The 
time  dependent  response  of  adhesively  bonded  joints  raises  questions 
regarding  their  long  term  reliability  under  creep  or  other  more  compli¬ 
cated  loading.  A  delayed  failure  (creep  rupture)  long  after  the 
initial  design  and  fabrication  process  is  possible.  Thus,  a  method  is 
needed  by  which  long  term  (a  few  years)  bonded  joint  response  can  be 
predicted  from  short  term  (a  few  days  at  most)  testing. 

Review  of  Literature 

The  eccentricity  of  the  applied  loads  as  well  as  the  bending  due 
to  the  nature  of  the  joint,  both  serve  to  complicate  the  stress  analysis 
of  the  single-lap  test  specimen.  The  first  and  probably  most  often 
quoted  paper  on  the  stress  analysis  of  the  single-lap  joint  was  pre¬ 
sented  by  Goland  and  Reissner  [34]  in  1944.  Their  solutions,  based  on 
the  bending  of  cylindrical  plates,  are  applicable  to  only  two  limiting 
cases:  (1)  when  the  adhesive  flexibility  can  be  neglected,  and  (2) 
when  the  adhesive  is  relatively  flexible.  The  solutions  are  based  on 
the  principle  of  minimum  potential  energy  and  are  restricted  to  linear 
isotropic  materials  and  identical  adherends.  Thermal  effects  were 
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neglected  in  the  solutions. 

The  stress  analysis  of  bonded  joints  has  received  a  great  deal  of 
attention  since  the  early  work  of  Go! and  and  Reissner.  Many  researchers 
have  attempted  to  predict  the  stresses  in  bonded  joints  using  closed 
form  analytical  solutions.  However,  the  resulting  equations  that  must 
be  solved  are  extremely  complicated  thus  requiring  that  simplifying  as¬ 
sumptions  be  made  as  in  the  Goland  and  Reissner  solutions.  Recently, 
the  finite  element  method  has  been  used  extensively  in  the  analysis  of 
bonded  joints  [18,19,35-37].  The  nonlinearity  of  the  adhesive  has  been 
taken  into  account  [18,19]  as  well  as  thermal  effects  [35].  The  time 
dependent  behavior  of  a  linear  viscoelastic  adhesive  has  also  been 
modeled  via  the  finite  element  method  [36,37]. 

The  common  denominator  between  all  of  the  bonded  joint  analyses 
performed  to  date,  whether  closed  form  or  finite  element,  is  that  they 
all  require  both  Young's  modulus  and  the  shear  modulus  for  the  adhesive 
layer  in  order  to  calculate  the  stress  distribution.  Young's  modulus 
is  usually  obtained  from  a  tensile  test  of  the  bulk  adhesive,  while  the 
shear  modulus  is  determined  via  a  single-lap  shear  test.  As  men¬ 
tioned  previously,  due  to  inherent  problems  with  the  single-lap 
specimen,  the  shear  modulus  may  be  more  accurate  if  obtained  from  a 
thick-adherend  short-overlap  shear  specimen. 

Efforts  have  been  made  to  determine  whether  the  mechanical 
properties  of  an  adhesive  in  a  joint  can  be  related  to  the  properties 
of  the  bulk  adhesive  [38-41].  Even  though  great  care  is  taken  to  cure 
the  adhesive  identically  for  both  situations,  the  thin  film  adhesive 
shear  properties  in  general,  cannot  be  related  to  the  bulk  tensile 


properties.  Although  this  is  the  case,  a  common  practice  in  the  past 
has  been  to  use  the  bulk  mechanical  properties  of  adhesives  for  the 
analysis  and  design  of  bonded  joints  [9]. 

Accurate  predictions  of  time  dependent  bonded  joint  response  have 
been  made  based  on  t!  time  dependent  bulk  tensile  adhesive  properties 
of  a  linear  viscoelastic  adhesive  [36,37].  However,  the  accuracy  of 
these  predictions  is  a  direct  consequence  of  the  linear  conditions  of 
small  stresses  and  strains.  Under  conditions  of  large  deformations,  as 
when  failure  is  approached,  this  technique  would  not  yield  satisfactory 
results.  In  the  large  strain  situation  both  the  bulk  and  the  thin  film 
material  properties  would  be  required  for  analysis. 

Objective  of  Overall  Research  Program 

The  present  investigation  represents  the  first  phase  of  a  multi¬ 
year  Virginia  Tech-NASA  Langley  collaborative  research  program  on  "The 
Viscoelastic  Characterization  and  Lifetime  Predictions  of  Structural 
Adhesives."  The  overall  objective  of  this  research  effort  is  the 
development  of  an  accelerated  characterization  procedure  which  would 
allow  the  prediction  of  failures  in  adhesively  bonded  joints  long  after 
the  initial  design  and  fabrication  process,  based  on  short  term  testing 
The  purpose  of  this  first  phase  was  to  study  the  uniaxial  tensile 
viscoelastic  behavior  of  two  commonly  used  structural  adhesives, 

FM-73*  and  FM-300*,  both  in  bulk  form.  A  nonlinear  viscoelastic 
characteri zation  procedure  was  to  be  applied  to  the  uniaxial  tensile 


♦Manufactured  by  the  American  Cyanamid  Company. 


creep  and  creep  recovery  data  for  both  materials.  Initially,  a  procedure 
due  to  Schapery  [29-32]  was  to  be  attempted  in  which  he  assumed  a 
creep  power  law.  However,  as  will  be  shown  subsequently,  the  power  law 
parameters  could  not  be  determined  for  each  of  our  particular  cases 
as  he  originally  suggested. 

The  Schapery  procedure  has  been  utilized  in  the  past  at  VPI  for 
adhesive  characterization.  Cartner  [25]  successfully  employed  the 
method  to  characterize  another  epoxy  structural  adhesive,  Metlbond 
1113-2,  in  bulk  form,  under  uniaxial  tensile  loading. 

Recent  research  by  Peretz  and  Weitsman  [33]  has  verified  the 
ability  of  the  Schapery  method  to  characterize  the  uniaxial  tensile 
creep  and  creep  recovery  response  of  FM-73,  in  bulk  form,  at  30°C 
(86°F). 

The  goal  of  the  current  research  effort  was  to  learn  how  to  pre¬ 
dict  the  long  term  uniaxial  tensile  creep  and  creep  recovery  response 
at  any  arbitrary  temperature  and  stress  level  from  short  term  data 
through  the  use  of  a  nonlinear  viscoelastic  mathematical  model  such  as 
that  of  Schapery  [29-32]  or  others. 

The  long  term  goal  would  be  of  course  to  utilize  the  time 
dependent  uniaxial  tensile  bulk  adhesive  properties  (determined  in  the 
present  investigation)  to  predict  time  dependent  bonded  joint  response 
under  linear  conditions.  This  could  be  done  using  a  procedure  similar 
to  that  employed  in  references  [36,37]. 

A  logical  extension  would  be  to  use  a  nonlinear  viscoelastic 
method  to  characterize  bonded  joint  response.  To  the  author’s  knowledge, 
a  nonlinear  viscoelastic  characterization  procedure  has  not  previously 


been  utilized  for  the  characterization  of  adhesively  bonded  joints. 

Upon  such  satisfactory  characterization  of  a  bonded  joint  the  analysis 
could  be  combined  with  a  time  dependent  failure  model  such  as 
developed  by  Zhurkov  [42]  or  Crochet  [43]  in  order  to  predict  long  term 
failures.  Such  a  procedure  would  give  a  methodology  with  which  the 
time  dependent  response  of  adhesively  bonded  joints  could  be  determined 
from  the  instant  of  initial  loading  until  failure  occurs,  thus 
accomplishing  the  overall  project  goal. 


Chapter  2 

APPROACHES  FOR  VISCOELASTIC  MODELING 


Linear  Viscoelasticity 

The  constitutive  equation  for  linear  elastic  materials  is  given 
by. 


eij  =  Sijklakl 


(1) 


where  S.^  are  the  compliances.  For  linear  viscoelastic  materials  un¬ 
der  creep  loading, 

£ij(t)  =  Si jkl (t)okl  (2) 

where  is  the  time  independent  stress  level  and  S^-j  are  the 
creep  compliances. 

For  more  general  loading  states,  one  may  express  the  strain 
response  as. 


eij(t^  =  Sijkl(t^°kl  H(t)  +  Sijkl(t  '  ll  Kl  H(t  '  tl 5  + 
which  may  be  generalized  by  the  following  Convolution  integral: 


£ 


d°kl(t> 

dx 


dx 


(4) 


When  the  material  is  unaffected  by  events  prior  to  t  =  0,  the  lower 
limit  of  the  integral  may  be  changed  to  zero. 

For  uniaxial  loading  and  no  prior  stress  history,  the  stress- 
strain  constitutive  equation  may  be  written  as, 


e(t) 


ft 

D(t  - 

.0 


x) 


dg(0 

dx 


dx 
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where  e(t)  is  the  time  varying  strain  output  for  an  arbitrary  time  vary¬ 
ing  stress  input,  a(t),  and  D(t)  is  the  creep  compliance.  Alternative¬ 
ly,  equation  (5)  may  be  written  as. 


e(t)  =  DQa(t)  +  f  AD(t  -  t)  — dx  (6) 

0  J  Q 

where  the  creep  compliance  has  been  separated  into  initial  and  transient 
components  such  that, 

D(t)  =  Dq  +  AD(t)  (7) 

The  Convolution  integral  of  equation  (4),  as  well  as  the  analogous 
equations  (5)  and  (6),  are  all  forms  of  the  Duhamel  or  Boltzmann 
superposition  integral. 

Time  Temperature  Superposition  Principle  (TTSP) 

The  aim  of  any  accelerated  characterization  scheme  is  that  one 
can  in  some  way  use  short  term  experimental  data  to  predict  long  term 
material  response.  The  Time  Temperature  Superposition  Principle  (TTSP) 
initially  proposed  by  Leaderman  [44]  has  been  widely  used  for  the 
accelerated  characterization  of  polymers.  The  basic  premise  of  the 
TTSP  is  that  compliance  curves  at  different  temperatures  are  of  the 
same  basic  shape,  but  only  shifted  in  time.  The  implication  of  the 
TTSP  is  that  compliance  data  at  several  different  temperatures  may  be 
shifted  horizontally  in  log  time  (vertical  shifting  may  also  be  re¬ 
quired)  to  produce  a  smooth  "master  curve"  approximating  the  compliance 
over  several  decades  of  time.  Although  it  has  be°'  su  :essfully 
employed  below  the  glass  transition  temperature,  T  ,  the  TTSP  is 
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rigorously  justifiable  only  above  the  T  .  It  should  be  noted  that  even 
though  it  is  not  a  requirement  of  the  theory,  the  TTSP  is  usually 
associated  with  linear  viscoelasticity  [45].  For  a  detailed  discussion 
of  the  TTSP  and  of  the  various  techniques  utilized  to  determine  the 
amount  of  horizontal  and  vertical  shifting  required  to  produce  a  smooth 
"master  curve"  the  reader  is  directed  to  the  work  of  Griffith  [45]. 

The  recent  research  of  Kenner,  Knauss,  and  Chai  [46]  employed  the 
TTSP  for  the  accelerated  characterization  of  FM-73.  Bulk  shear  creep 
compliance  curves  for  several  temperatures  (above  and  below  the  T  ) 
were  shifted  both  vertically  and  horizontally  to  form  a  "master  curve" 
at  20.5°C,  extending  over  sixteen  decades  of  log  time.  The  applied 
torsional  stresses  were  sufficiently  small  that  linear  viscoelastic 
behavior  was  assumed. 

Time  Stress  Superposition  Principle  (TSSP) 

Analogous  to  the  TTSP  is  the  Time  Stress  Superposition  Principle 
(TSSP)  which  may  be  given  by  the  basic  equation, 

D(a,T)  =  D0(o)  +  baA0(c)  (8) 

where  ba  is  the  vertical  shift  factor,  ^  is  the  reduced  time  given 

by  and  a.  is  the  horizontal  shift  factor  due  to  the  stress  level. 

aa  J 

The  implication  of  the  TSSP  is  that  isothermal  compliance  data  at 
various  stresses  (for  nonlinearly  viscoelastic  materials)  may  be 
shifted  both  horizontally  and  vertically  to  form  a  "master  curve"  pre¬ 
dicting  long  term  compliance  based  on  short  term  testing. 
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Time  Temperature  Stress  Superposition  Principle  (TTSSP) 

The  similarities  between  the  TTSP  and  the  TSSP  are  apparent.  In 
the  TTSP,  temperature  is  the  accelerating  parameter  whereas  in  the  TSSP 
stress  is  the  accelerating  parameter.  It  would  seem  logical  that  an 
accelerated  characterization  technique  which  would  simultaneously 
employ  both  stress  and  temperature  as  accelerating  parameters  would  be 
possible.  Indeed,  elements  of  the  TTSP  and  the  TSSP  have  been  utilized 
to  form  the  Time  Temperature  Stress  Superposition  Principle  (TTSSP). 
Again  the  reader  is  directed  to  the  work  of  Griffith  [45],  which  dis¬ 
cussed  the  determination  of  the  appropriate  horizontal  and  vertical 
shift  functions  for  the  TTSSP. 

The  application  of  the  TTSP  as  well  as  the  TSSP  and  the  TTSSP  is 
graphical  in  nature  requiring  a  large  data  base  and  thereby  creating 
a  tedious  and  time  consuming  approach. 

Analytical  Nonlinear  Viscoelastic  Approaches 

Several  techniques  have  been  developed  to  account  for  nonlinear 
viscoelatic  behavior  [29-32,47-60].  The  theories  of  Green  and  Rivlin 
[47],  Green,  Rivlin,  and  Spencer  [48],  and  Green  and  Naghdi  [49],  are 
all  triple  integral  representations  which  are  very  difficult  to  solve. 
Also  difficult  to  solve  are  the  coupled  first  order  nonlinear  dif¬ 
ferential  equations  of  Walker  [50],  and  Krempl  [51-56].  Two  relatively 
simple  nonlinear  theories  are  those  of  Schapery  [29-32],  and  Findley 
[57-60].  Although  essentially  an  empirical  curve  fitting  procedure, 
the  Findley  technique  is  simple  to  utilize  and  has  been  successfully 
employed  for  long  term  creep  predictions. 


Because  elements  of  both  the  Schapery  and  Findley  methods  will 
be  used  subsequently  to  analyze  data,  these  two  approaches  will  be 
briefly  reviewed  here. 


Schapery  Method 

The  first  nonlinear  viscoelastic  characterization  procedure  em¬ 
ployed  in  the  current  work  was  developed  by  R.  A.  Schapery  [29,30]  in 
the  late  1960’s.  The  derivation  of  the  theory  is  too  lengthy  and  de¬ 
tailed  to  discuss  here,  but  is  reviewed  and  discussed  in  detail  by 
Hiel  [61], 

Schapery1 s  nonlinear  viscoelastic  stress-strain  equation  for  one 
dimensional  loading  may  be  written  as  [30], 


e(t)  *  90D0°U)  +  g1 


ft 

AD (v  -  Ip' ) 

•*0 


b  [g2a(T)]dT 


(9) 


where  0Q  and  aD(*)  are  the  initial  and  transient  components  of  the 
linear  viscoelastic  creep  compliance.  The  reduced  time  parameters. 


v  =  i'(t)  = 


<P'  =  V  (t) 


(10) 


are  a  function  of  the  stress  dependent  shift  function,  a  .  The  quanti¬ 
ties,  g  ,  g-j  and  g^,  are  stress  dependent  properties  that  represent  the 
nonlinear  nature  of  a  material  and  must  be  determined  experimentally. 
When  gQ  =  g-]  =  a^  =  1 ,  equation  (9)  reduces  to  the  Modified  Superposi¬ 
tion  Principle  (MSP)  proposed  by  Leaderman  [44],  When  all  stress 
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dependent  properties,  g  ,  g-j ,  g2  and  ag ,  are  taken  to  be  equal  to  one, 
equation  (9)  reduces  to  the  Boltzmann  superposition  integral  of  equa¬ 
tion  (6).  Though  not  specifically  stated  by  Schapery,  his  approach 
represented  by  equation  (9)  is  a  mathematical  statement  of  the  Time 
Stress  Superposition  Principle  (TSSP),  which  was  discussed  previously. 

All  of  the  material  properties  in  equation  (9)  can  be  evaluated 
from  the  data  obtained  from  a  creep  and  creep  recovery  test  [32].  In 
this  test,  a  specimen  is  subjected  to  a  constant  stress,  a  ,  which  is 
maintained  until  time  t-j ,  then  removed  (Figure  1).  The  stress  input 
for  a  uniaxial  creep  and  creep  recovery  test  is  given  by, 

a(t)  =  a0H(t)  -  aQ H(t  -  t})  (11) 

where  aQ  is  the  time  independent  input  stress  level,  and  H(t)  is  the 
Heaviside  step  function  defined  as, 


H(t)  =  0,  when  t  <  0 
H(t)  =  1 ,  when  t  >  0 


(12) 


Substitution  of  equation  (11)  into  equation  (9)  can  be  shown  to 

yield, 

ec(t)  =  [gQ0o  +  g1g2AD(|— )]  ctq  (13) 

a 

and 

t, 

er(t)  =  [g2AD(^-  +  t  -  tj)  -  g24D(t  -  tj)]  aQ  (14) 

“  o 

for  the  creep  and  creep  recovery  strain,  respectively. 

If  it  is  assumed  that  the  compliance  follows  a  power  law  in  time, 
as  do  many  polymers,  the  transient  portion  may  be  written  as, 

aDU)  =  C*n 


(15) 


If  equation  (15)  is  substituted  into  equations  (13)  and  (14),  the 
expressions  for  creep  and  creep  recovery  strain  become, 

ec(t)  =  goDo°o  +  Caogl92(r)  (16) 

a 

and 

er(t)  ■  p-  [(1  +  a0x)n  -  (a0x)n]  (17) 

where,  C  and  n  are  material  properties  which  are  invariant  with  stress, 

t  -  t, 

X  =  — r— 1  (18) 

is  a  nondimensional  time,  and 

Ae1  =  e ( 1 1  )  -  eQ  (19) 

In  Equation  (19),  Ae^  represents  the  transient  component  of  the  creep 
strain  just  prior  to  the  removal  of  the  stress  at  t  =  t-j . 

Schapery's  approach  was  derived  from  the  theory  of  irreversible 
thermodynamics.  As  a  result,  each  of  the  stress  dependent  properties 
in  equation  (9)  has  a  thermodynamic  origin.  Schapery  [30]  has  shown 
that  changes  in,  gQ,  g^ ,  and  g^,  reflect  dependence  of  the  Gibbs  free 
energy  on  the  applied  stress,  while  changes  in  a  arise  due  to  a 

3 

similar  stress  dependence  of  both  entropy  production  and  free  energy. 

The  greatest  advantage  of  the  Schapery  approach  over  the  nonlinear 
theories  mentioned  previously  (excluding  Findley)  is  the  relative  sim¬ 
plicity  of  the  calculations.  This  is  a  result  of  the  simple,  single 
integral  form  of  equation  (9).  An  additional  attribute  of  the  method 
is  that,  once  the  nonlinear  parameters  are  determined  from  tensile 
creep  and  creep  recovery  using  equations  (16)  and  (17),  they  can  be 
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utilized  in  the  more  general  equation  (9)  to  model  a  wide  variety  of 
loading  and  unloading  situations. 


Findley  Method 

A  nonlinear  viscoelastic  characterization  method  studied  exten¬ 
sively  by  Findley  [57-60]  was  also  incorporated  into  the  current 
analysis.  The  basic  concept  behind  the  Findley  analysis  is  that  for 
any  given  creep  load,  the  specimen  strain  is  given  by, 

e ( t )  *  eQ  +  m  tn  (20) 

where  eQ,  m,  and  n  are  material  properties.  Further,  the  assumptions 
are  made  that 

n  *  constant,  independent  of  stress  level 


,  *  e'  sinh  a/a 

oo  e 

(21) 

m  *  m'  sinh  a/a 

m 

(22) 

where  e',  a  ,  m' ,  and  a  are  material  constants  for  any  given  tempera- 
o  e  m 

ture,  moisture  level,  etc.  The  nonlinear  effect  of  stress  is  accounted 
for  by  the  hyperbolic  sine  terms. 

It  has  been  shown  by  Schapery,  that  his  expression  for  creep 
strain  given  by  equation  (16)  is  equivalent  to  the  Findley  creep  ex¬ 
pression  of  equation  (20)  if 


and 


sinh  a/a 


(23) 


g1g2  ^  sinh  a/am 


(24) 


Determination  of  the  Schapery  Parameters 


If  a  viscoelastic  material  is  to  be  characterized  by  equations 
(16)  and  (17),  seven  material  parameters  must  be  evaluated  (n,  C,  DQ, 
g  ,  g-j ,  gg.  and  a^).  Schapery  suggested  a  graphical  procedure  to 
determine  these  seven  parameters  [32].  A  brief  description  of  his 
procedure  follows. 

First,  the  exponent,  n,  is  determined  from  linear  creep  recovery 
data.  Taking  the  logarithm  of  equation  (17)  gives, 

Ae, 

log  er  =  log  er  -  log  —  *  log  [(1  +  aax)n  -  (a^x)"]  (25) 

where  log  er  is  defined  as  the  shifted  recovery  strain.  When  g-j  =  aa  = 
1  (linear  range),  a  family  of  double-logarithmic  curves  of  [(1  +  x)n 
-  x"]  versus  X  can  be  plotted  for  several  values  of  n  (usually 
0  <  n  <  0.5).  Experimentally  known  values  of  log  er  can  be  plotted 
against  X  in  the  linear  range  and  the  shape  of  the  resulting  curve  can 
be  matched  to  one  of  the  aforementioned  family  of  curves,  thereby 
determining  n.  Next,  log  for  each  stress  level  in  the  nonlinear 
range  can  be  put  on  the  same  graph  with  the  linear  data.  Since  n  is 
independent  of  stress  level,  all  the  curves  of  log  cr  for  each  stress 
level  can  be  shifted  to  form  a  single  continuous  creep  recovery  "master 
curve."  The  necessary  horizontal  shifting  for  each  stress  level  is 


log  a^  while  the  required  vertical  shifting  for  each  stress  level 


represents  log 


£fl 

9l 


Let  equation  (16)  be  represented  by, 
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where  *  g  D  a  is  the  initial  strain  and  the  creep  coefficient,  C', 
0  0  0  0 


is  given  by, 

.  C°p3l92 


(27) 


The  quantities  eQ  and  C'  can  be  obtained  for  each  stress  level  by 
solving  equation  (26)  at  any  two  times,  tfl  and  tfa. 

At  this  point  the  transient  creep  strain,  Ae-j ,  can  be  determined 
for  each  stress  level  as, 


Ae1-C'(t1)n  (28) 

Ae, 

As  the  quantity  is  known,  g-j  may  now  be  easily  determined. 

For  linear  data  equation  (27)  reduces  to,  C1  =  CaQ,  thus  the 
parameter  C  may  be  evaluated.  The  parameter  g2  can  be  determined  for 
all  nonlinear  stress  levels  from  equation  (27). 


The  remaining  parameter,  gQ,  is  evaluated  for  each  stress  level 


9°  Da 
0  0 


The  above  graphical  procedure  is  tedious,  time  consuming,  and 
quite  subjective.  For  this  reason,  a  least  squares  computer  based 
procedure  has  been  developed  at  VPI  by  Hi  el  [61]  and  Bertolotti  [62]. 
The  method  uses  a  subroutine  labeled  "ZXSSQ"  developed  by  IMSL,  Inc. 


for  an  IBM/SINGLE  computer.  The  details  of  this  program  are  discussed 


in  the  following  chapter. 


Chapter  3 

APPROACHES  TO  DETERMINING  MATERIAL  PARAMETERS 


Computer  Modified  Schapery  Method 

Recall  that  the  expressions  for  uniaxial  tensile  creep  and  creep 
recovery  strain  were  given  by, 

£c<4)  •  9oVo  +  Coo9l92(r')"  (30) 

a 

and 

c_(t)  •  ^-£(1  +  a  xf  -  (a  x)n]  (31) 

”9]  a  a 

In  order  to  find  the  seven  material  parameters  in  equations  (30)  and 
(31)  Schapery  [32]  suggested  that  creep  and  creep  recovery  tests  in 
both  the  linear  and  nonlinear  range  of  the  material  were  required.  A 
typical  testing  program  should  consist  of  at  least  one  test  in  the 
linear  range  and  several  tests  in  the  nonlinear  range.  The  above 
series  of  tests  would  be  performed  under  isothermal  (constant  tempera¬ 
ture)  conditions.  It  should  be  emphasized  that  in  this  procedure  the 
power  exponent  is  determined  from  creep  recovery  data.  The  reason  for 
this  is  that  transient  strains  are  only  present  in  recovery  data  and 
these,  of  course,  correspond  to  the  second  term  in  equation  (30). 

The  computer  based  method  of  determining  the  unknown  parameters 
of  equations  (30)  and  (31)  developed  by  Heil  and  Bertolotti  [62]  in¬ 
corporates  a  very  powerful  subroutine  labeled  "ZXSSQ"  developed  by 
IMSL,  Incorporated.  This  subroutine  uses  a  finite  difference  Levenberg 
Marquardt  algorithm  to  find  the  minimum  sum  of  squares  of  the  dif¬ 
ference  between  a  number  of  experimental  data  points  and  a  user  defined 
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function.  Convergence  is  satisfied  and  execution  is  terminated  if  on 
two  successive  iterations,  the  parameter  estimates  agree,  component  by 
component,  to  three  digits. 

This  subroutine  allows  the  user  to  define  a  curve  that  is  to 
model  the  data.  The  unknown  parameters  are  evaluated  in  such  a  way 
as  to  minimize  the  error  between  each  experimental  data  point  and  the 
defined  curve.  A  more  detailed  discussion  of  this  subroutine  is  given 
later  in  this  chapter. 

The  computer  program  was  set  up  to  analyze  one  set  of  creep  and 
creep  recovery  data  for  each  run.  The  first  computer  run  is  made 
using  the  linear  creep  and  creep  recovery  data.  In  the  linear  range, 
by  virtue  of  the  fact  that  g-j  =  aa  =  1 ,  the  creep  recovery  equation 
(31)  reduces  to, 

Cr(t)  =  Ae,[(l  +  X)n  -  Xn]  (32) 

Subroutine  "ZXSSQ"  is  then  used  to  determine  the  best  fit  for  ^  and  n. 
Letting  gQ  =  g-|  =  g£  =  =  1  in  equation  (30),  we  get, 

ec(t)  =  Vo  +  Cao(t)"  (33) 

where  n  is  now  known.  Equation  (33)  can  be  written  as, 

€c(t)  =  +  C't"  (34) 

Letting  subroutine  "ZXSSQ"  operate  on  equation  (34),  yields  best  fit 

values  for  e  and  C' ,  where 
o 

t  -  Da 
o  oo 

and 

c  ■  V 


(36) 


Thus  from  analysis  of  the  linear  data  the  three  parameters,  DQ,  C  and 
n  are  now  known.  These  three  parameters  are  assumed  to  be  invariant 
with  stress,  and  thus  will  remain  constants  throughout  the  nonlinear 
analysis  to  follow. 

The  computer  program  is  now  executed  with  the  first  set  of  non¬ 
linear  creep  and  creep  recovery  data.  Subroutine  "ZXSSQ"  operates  on 

equation  (31)  with  n  as  a  known  quantity  and  determines  the  quantities 

Ae, 

-jj-  and  a^.  Subroutine  "ZXSSQ"  then  operates  on  equation  (34)  with  n 
as  a  known  quantity  and  determines  the  best  fit  for  and  C',  where 

£  =  5  Do  (37 

0  0  0  0 


Realizing  that. 


C°o9192 


Ac,  =  C*  (t,  )n 


g-,  can  be  easily  determined  by, 


«i  ■  i4f]  •  c'(t')n 

From  equation  (38),  g^  can  be  determined  as, 

C'U/ 

32  '  Cog, 


The  last  parameter  g  can  be  determined  from  equation  (37), 
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The  computer  program  is  then  reexecuted  for  each  set  of  nonlinear 
creep  and  creep  recovery  data. 

Illustration  of  Subroutine  "ZXSSQ" 

The  recovery  equation  (32)  will  be  used  as  an  example  to  illustrate 
the  numerical  process  involved  in  evaluating  the  unknowns  and  n. 
Subroutine  "ZXSSQ"  requires  that  the  right  hand  side  of  equation  (32) 
be  written  as, 

X(l)*((l  ,+XLAB( I )**X(2)-XLAB( I )**X(2) )  (43) 

Let, 

F(I)=REC(I)-X(1  )*( (1  .+XLAB(I )**X(2)-XLAB(I )**X(2) )  (44) 

where, 

X(1 )  and  X(2)  —  parameters  to  be  evaluated  (X(l)  replaces 
Ae-j  and  X(2)  replaces  n) 

REC(I)  --  array  containing  the  experimental  creep  recovery 

strai ns 

XLAB(I)  --  array  containing  the  reduced  time  values 

F(I)  —  array  containing  the  error  between  the  experi¬ 

mental  creep  recovery  strains  and  the  recovery 
equation  (43) 

A  DO  LOOP  around  equation  (44)  evaluates  the  values  of  the  array, 

F(I).  Each  of  these  values  is  squared  and  all  elements  are  summed. 

Thus  for  every  value  of  X(l)  and  X(2)  there  is  an  associated  value  of 

2 

the  sum  of  squres  error,  zR  .  If  a  series  of  values  of  X(l)  and  X ( 2 ) 

2 

and  their  associated  error,  zR  ,  are  plotted  on  a  3-D  graph,  an  error 
surface  is  formed  (Figure  2). 


The  unknown  parameters  X(.l)  and  X(2)  must  be  given  initial  values 


by  the  user.  Subroutine  "ZXSSQ"  begins  searching  for  the  best  fit 
values  of  X(l)  and  X(2)  by  examining  the  initial  guesses.  The  sub¬ 
routine  calculates  the  value  of  the  normal  to  the  gradient  for  the 
initial  values  of  X(l)  and  X(2),  and  finds  the  direction  of  steepest 
descent.  New  values  of  X(l)  and  X(2)  are  found  in  the  direction  of  the 
steepest  descent  line,  and  the  normal  to  the  gradient  is  calculated 
again.  This  process  continues  until  the  parameter  estimates  of  two 
successive  iterations  agree  to  three  digits.  The  procedure  just 
described  is  illustrated  in  Figure  2. 

Computer  Modified  Findley  Approach 

As  will  be  discussed  in  a  later  chapter,  the  computer  method 
described  above  did  not  always  result  in  stable  values  for  n  and  C. 

As  a  result,  a  Findley  based  approach  to  the  determination  of  these 
parameters  which  uses  only  creep  data  was  utilized.  The  computer 
program  to  so  determine  n  and  C  was  developed  by  Yen  [63]  using  the 
same  subroutine,  "ZXSSQ."  This  modified  approach  is  in  fact  equivalent 
to  that  employed  by  Peretz  and  ’^itsman  [33]  and  determines  the  three 
parameters,  C,  n  and  D  ,  from  linear  creep  data  as  discussed  below. 

Recall  that  the  creep  strain  can  be  written  as, 

ec(t)  =  eQ  +  C'tn  (45) 

Letting  subroutine  "ZXSSQ"  perform  a  three  parameter  fit  on  equation 
(45),  will  yield  the  best  fit  values  of  e  ,  C'  and  n,  where 
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L\*. 


and 

C'  =  CaQ  (47) 

Thus  from  analysis  of  the  linear  creep  data  only,  the  three  parameters, 
Dq,  C  and  n  can  be  determined.  These  parameters  then  remain  invariant 
throughout  the  nonlinear  analysis. 

The  nonlinear  analysis  remains  as  was  previously  detailed. 

A£] 

Specifically,  the  quantities  — —  and  a  are  determined  from  a  two 

g-|  a 

parameter  fit  of  equation  (31)  and  the  quantities  zQ  and  C'  are  fit 
from  equation  (34).  The  three  nonlinear  parameters,  g  ,  and  are 
then  determined  from  equations  (42),  (40)  and  (41),  respectively. 

The  inconsistencies  between  the  two  procedures  and  the  results 
obtained  will  be  discussed  in  detail  in  Chapter  5. 


Chapter  4 

EXPERIMENTAL  PROCEDURES 

Cure  Cycle  and  Specimen  Preparation 

The  two  adhesives  investigated  in  this  study  were  FM-73  and 
FM-300.  Both  adhesives  are  produced  by  the  American  Cyanamid  Company 
of  Havre  de  Grace,  Maryland.  Both  FM-73  and  FM-300  are  manufactured 
and  distributed  in  the  form  of  a  thin  film  supported  by  a  polyester 
carrier.  These  thin  films  (plies)  of  adhesive  material  were  formed 
into  bulk  adhesive  panels  at  NASA-Langley  Research  Center.  The  cure 
cycles  employed  by  NASA-Langley  are  given  in  Table  1. 

The  final  bulk  adhesive  panels  were  approximately  12  x  12  x  0.05 
inches.  Dog  bone  specimens  whose  dimensions  are  given  in  Figure  3, 
were  machined  from  each  panel. 

Equi pment 

All  tests  were  performed  on  an  ATS  (Applied  Test  Systems)  Model 
2330  lever  arm  creep  test  machine  with  automatic  draw  head  and 
releveler.  Loading  or  unloading  required  approximately  10  seconds. 
Since  this  time  was  less  than  IX  of  the  total  test  time,  the  experiment 
was  considered  to  be  a  good  approximation  of  instantaneous  (step) 
loading. 

The  creep  machine  was  equipped  with  an  ATS  series  2912  oven  and 
a  series  230  temperature  controller.  This  control  unit  maintained  an 
oven  temperature  within  t2°F  (±1°C)  of  the  desired  temperature. 


Table  1.  Adhesive  Cure  Cycles 


FM-73 

7  plies  --  10  mils  per  ply 
full  vacuum 


Temperature 
1 50°F  (65.5°C) 
200°F  (93.3°C) 
250°F  (121. 1°C) 


Time 


1  hour  at  each 
temperature 


Vi  ton  Dam  &  Plus  —  Nylon  Sheet  each  side 
Sealant  on  base  to  retain  resin  (adhesive) 


FM-300 


6  plies  --  16  mils  per  ply 
full  vacuum 

Temperature  Time 

250°F  (121. 1°C) 

300°F  (148. 9°C) 

350°F  (176. 7°C) 


1  hour  at  each 
temperature 


Viton  Dam  &  Plus  -  Nylon  Sheet  on  each  side 


igure  3.  Tensi 


Temperature  was  monitored  by  thermocouples  and  a  Doric  41 2A  Trendicator 
digital  thermometer. 

The  more  desirable  procedure  of  using  back  to  back  gages  was  not 
utilized  due  to  the  small  thickness  (and  the  corresponding  cross 
sectional  area)  of  the  test  specimens.  It  was  felt  that  the  back  to 
back  gage  arrangement  would  produce  a  significant  stiffening  of  the 
specimen.  The  small  specimen  thickness  also  tends  to  minimize  the  in¬ 
duced  bending  effects  caused  by  eccentricities  in  loading.  It  is  this 
induced  bending  which  usually  necessitates  the  use  of  the  back  to  back 
gage  arrangement.  Caplan  [64]  conducted  tests  on  polycarbonate  using 
the  back  to  back  gage  arrangement.  He  found  the  difference  between  the 
two  strain  gage  readings  to  be  negligible.  It  was  thus  determined  that 
back  to  back  gages  were  both  unnecessary  and  undesirable  in  the  present 
situation.  For  all  tests,  strain  was  measured  using  a  single  strain 
gage.  These  gages  were  Micro-Measurements  EP-08-1 25BB-1 20 ,  120  ohm 
gages  with  pre-attached  lead  wires.  Gages  with  pre-attached  lead  wires 
were  used  to  eliminate  soldering  directly  on  the  specimen  surface. 
Surface  preparation  was  performed  according  to  recommended  Micro- 
Measurements  procedures.  Gages  were  bonded  with  M-Bond  600  adhesive 
and  cured  according  to  the  manufacturer's  specifications.  An  un¬ 
strained  specimen  with  a  "dummy"  (compensating)  gage  was  used  in  a  half¬ 
bridge  arrangement  to  compensate  for  strain  due  to  fluctuations  in 
temperature.  A  voltage  of  2V  was  used  to  minimize  induced  strains 
due  to  gage  heating  effects  in  the  120  ohm  gages.  A  2120  Vishay  System 
was  used  to  condition  the  strain  gage  output.  The  conditioned  strain 
data  was  read  from  an  MTS  408.31  digital  voltmeter. 


The  specimens  were  pin  loaded  using  serrated  metal  grips  bolted 


around  the  specimen  ends. 

Moisture  Content 

Moisture  content  is  known  to  have  a  significant  effect  on  the 
strength  and  stiffness  of  polymers.  Absorption  of  moisture  is  known  to 
have  an  effect  similar  to  increasing  the  temperature  [65].  It  has  been 
shown  [32]  that  even  a  small  change  in  humidity  can  affect  specimen 
response  appreciably.  Therefore,  moisture  content  could  be  a  major 
factor  in  creep  studies.  However,  for  the  purposes  of  the  present 
study,  it  seemed  reasonable  to  characterize  the  bulk  adhesives  in  the 
as-received  condition.  Therefore  all  specimens  were  stored  and  tested 
under  normal  laboratory  conditions. 

Mechanical  Conditioning 

It  became  apparent  very  early  in  the  study  that  repeatable  creep 
and  creep  recovery  curves  could  not  be  obtained  from  a  specimen  which 
had  not  been  "conditioned."  It  was  also  observed  that  a  specimen  which 
had  not  been  conditioned  would  not  return  to  a  state  of  zero  strain 
after  the  load  was  removed.  The  reason  for  this  permanent  strain 
(damage)  was  assumed  to  be  the  occurrence  of  crack  and  damage  growth 
during  the  loading  phase  [32,33].  It  was  felt  that  if  the  specimen 
was  repeatedly  loaded  and  unloaded  (mechanically  conditioned),  a 
constant  damage  state  would  be  reached  and  repeatable  results  could 
then  be  obtained  [32,33],  This  was  found  to  be  the  case. 

Following  a  procedure  similar  to  Peretz  and  Weitsman  [33],  the 
gaged  dog  bone  specimens  were  mechanically  conditioned  by  subjecting 


them  to  several  loading  and  unloading  cycles  at  approximately  80%  of 
their  ultimate  strength.  A  particular  creep  load  was  applied  to  a 
specimen  and  maintained  for  approximately  20  seconds  and  then  removed. 
This  process  was  repeated  about  ten  times  or  until  satisfactory  stress- 
strain  reproducibility  was  attained.  It  should  be  noted  that  even 
though  a  specimen  may  have  been  conditioned  previously,  if  the  specimen 
was  to  be  used  at  a  higher  temperature  than  that  at  which  it  was 
conditioned,  the  conditioning  process  was  repeated.  This  was  done 
because  polymeric  materials  often  "heal"  themselves  upon  unloading  and 
the  healing  process  is  accelerated  by  an  increase  in  temperature. 

Testing  Program 

A  series  of  creep  and  creep  recovery  tests  were  performed  at 
various  temperatures  and  stress  levels.  The  test  temperatures  selected 
were,  room  temperature,  40°C  (104°F),  55°C  (1 31 °F) ,  70°C  (158°F), 

85°C  (185°F),  and  1Q0°C  (212°F).  Similar  to  Peretz  and  Weitsman  [33], 
the  lowest  stress  level  used  at  each  temperature  was  assumed  to  be  in 
the  linear  range  of  the  material,  even  though  FM-73  and  FM-300  are 
linear  only  for  exceedingly  small  stresses  and  strains,  if  they  are 
linear  at  all.  The  assumption  of  linearity  was  necessary  because  of 
the  difficulty  of  obtaining  accurate  strain  values  associated  with 
small  stress  levels. 

Although  a  multitude  of  creep  and  creep  recovery  tests  have  been 


performed  to  date,  and  the  results  have  been  reported  and  discussed 
by  many  authors,  there  is  no  consensus  at  the  present  time  as  to  the 
ideal  length  of  a  creep  and  creep  recovery  experiment.  Lou  and 


Schapery  £32]  used  a  test  program  consisting  of  one  hour  of  creep  and 
two  hours  of  recovery.  Peretz  and  Weitsman  [33],  with  the  aid  of  a 
data  acquisition  system  were  able  to  perform  the  nonlinear  viscoelastic 
characterization  of  FM-73  by  collecting  only  15  minutes  of  creep  data 
and  5  minutes  of  recovery  data.  Based  upon  creep  tests  performed  by 
Cap! an  [64]  on  polycarbonate,  it  was  decided  that  a  test  program  con¬ 
sisting  of  30  minutes  of  creep  and  60  minutes  of  creep  recovery  would 
be  sufficient  for  the  present  study. 

Inasmuch  as  the  purpose  of  the  present  effort  was  to  determine 
the  feasibility  of  using  a  nonlinear  viscoelastic  characterization 
procedure  for  FM-73  and  FM-300  and  not  to  determine  statistically  valid 
properties,  the  testing  program  needed  to  determine  the  nonlinear 
parameters  described  in  the  preceding  chapters  was  conducted  on  a 
single  specimen.  However,  the  two  short  term  experiments  performed 
to  predict  the  creep  and  creep  recovery  response  of  FM-300  under  an 
arbitrary  temperature  and  stress  situation  were  performed  with  a  "fresh 
specimen,  as  were  the  two  long  term  creep  experiments. 


Chapter  5 

RESULTS  AND  DISCUSSION 


FM-73  Adhesive  Characterization 

The  modified  epoxy  adhesive,  FM-73,  exhibited  viscoelastic  creep 
and  creep  recovery  phenomena  at  all  temperature  levels.  Isochronous 
stress-strain  results  for  30°C  (86°F)  are  shown  in  Figure  4.  The 
lowest  stress  level  employed,  493  psi  (3.45  MPa),  was  assumed  to  be 
the  upper  limit  of  linear  viscoelastic  material  response.  The  reason 
for  this  assumption  of  linearity  was  explained  in  the  preceding  chapter 
Since  this  initial  stress  level  is  less  than  10%  of  the  ultimate  stress 
of  the  material  the  assumption  seems  justifiable. 

The  Schapery  parameters  necessary  to  characterize  the  visco¬ 
elastic  response  of  FM-73  at  30°C  (86°F)  were  obtained  via  the  computer 
approach  described  in  Chapter  3.  The  stress  invariant  parameters,  n, 

C,  and  Dq,  were  determined  to  be,  0.151,  0.219,  and  1.569,  respec¬ 
tively.  The  stress  dependent  material  parameters,  gQ,  g-|  ,  g2»  and  a^, 
are  plotted  in  Figures  5-7.  With  the  seven  Schapery  parameters  known, 
the  complete  creep  and  creep  recovery  response  of  FM-73  at  30°C  (86°F) 
can  be  described  by  equations  (16)  and  (17)  respectively. 

The  vertical  shift  parameter,  (Ae-j/g^),  is  plotted  as  a  function 
of  stress  in  Figure  8.  Use  of  this  parameter  greatly  facilitates  the 
application  of  equation  (17)  to  describe  the  recovery  response.  As 
mentioned  previously,  the  creep  response  given  by  equation  (16)  can  be 
stated  in  a  simpler  fashion  by  equation  (26).  The  instantaneous  creep 
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strain,  eQ ,  and  the  creep  coefficient,  C',  both  plotted  as  a  function 
of  stress  level,  are  given  in  Figures  9  and  10,  respectively. 

A  comparison  of  experimental  creep  data  and  the  corresponding 
Schapery  fit  of  the  data  for  several  stress  levels  at  30°C  (86°F)  is 
given  in  Figure  11.  A  similar  comparison  for  recovery  data  is  given  in 
Figure  12.  Inspection  of  Figure  12  shows  that  the  Schapery  fit  yields 
strain  values  lower  than  the  experimental  values  as  the  length  of 
recovery  time  increases.  This  trend  in  the  Schapery  recovery  curve 
fit  was  apparent  throughout  the  current  research. 

A  master  creep  recovery  curve  for  FM-73  at  30°C  (86°F)  is  given  in 
Figure  13.  This  "master  curve"  is  included  here  to  illustrate  the  ac¬ 
curacy  with  which  the  computer  based  procedure  determines  the  hori¬ 
zontal  and  vertical  shift  factors,  a^  and  (ac-j /g^ )  -  In  the  past,  when 
the  Schapery  procedure  was  performed  graphically,  the  shift  factors,  a 
and  ( a e i / 9 1 ) ,  were  defined  as  the  amount  of  horizontal  and  vertical 
shifting,  respectively,  necessary  to  form  a  smooth  master  curve  from 
the  log-log  plot  of  recovery  strain  vs.  x  for  the  various  stress  levels. 
The  master  curve  of  Figure  13  was  generated  in  the  reverse  manner.  The 
computer  determined  values  for  a^  and  ( a e ^ /g i )  were  used  to  shift  the 
recovery  strain  vs.  a  curves  for  the  various  stress  levels.  Insnection 
of  Figure  13  reveals  that  this  shifting  did  result  in  a  smooth  creep 
recovery  master  curve,  thus  illustrating  the  accuracy  of  the  computer 
based  technique. 

It  was  attempted  at  this  point  to  perform  the  same  characteri za- 
tion  procedure  as  described  above,  for  FM-73  at  40°C  (104°F).  A 
complete  set  of  creep  and  creep  recovery  data  was  collected  at  this 
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Figure  11.  Experimental  creep  data  (symbols)  and  the  correspond 
ing  Schapery  fit  of  the  data  (solid  lines)  for 
several  stress  levels  at  T  =  30 0 C  (86’F). 


temperature.  The  data  was  analyzed  and  the  seven  material  parameters 
were  evaluated.  Plots  of  these  parameters  obtained  at  this  temperature 
exhibited  large  amounts  of  scatter,  so  much  scatter  in  fact  that 
fitting  a  curve  to  the  results  was  not  possible.  Upon  further  investi¬ 
gation  it  was  discovered  that  the  creep  and  creep  recovery  strain 
values  for  any  particular  stress  level  at  this  temperature  could  not  be 
reproduced  with  the  desired  accuracy. 

There  are  many  possible  reasons  that  could  account  for  the 
scatter  mentioned  above.  These  might  be  enumerated  as; 

•  Fluctuations  in  temperature  from  test  to  test  or  during  a  test 

•  Absorption  or  evaporation  of  moisture  between  tests  or  during 
a  test 

•  Experimental  error 

•  Data  reduction  error 

In  fact,  the  scatter  might  be  attributable  to  a  combination  of  all  of 
these  items.  For  example,  it  is  well  known  that  moisture  may  lower 
the  Tg  by  a  substantial  amount.  If  such  were  the  case,  test  tempera¬ 
tures  of  50  or  60°C  could  be  in  or  near  the  transition  region  of  the 
material  where  the  modulus  varies  Dy  a  large  amount,  giving  rise  to 
errors  in  measurement.  Also,  as  a  test  progresses,  the  specimen  is 
continually  drying  out,  thereby  essentially  raisinq  the  transition 
point  and  inducing  further  error.  In  addition,  the  Schapery  method 
requires  high  accuracy  in  measuring  transient  response.  Thus,  the 
problems  might  be  attributable  to  the  accuracy  involved  in  obtaining 
recovery  strains  and  the  data  reduction  related  thereto. 


The  modified  epoxy  adhesive,  FM-300,  was  found  to  exhibit  visco¬ 
elastic  creep  and  creep  recovery  phenomena  at  all  temperature  levels. 
Isochronous  stress-strain  results  for  FM-300  at  26°C  (79°F)  are  given 
in  Figure  14.  For  reasons  presented  previously,  the  lowest  stress 
level  employed,  452  psi  (3.17  MPa),  was  assumed  to  be  the  upper  limit 
of  linear  response  for  each  of  the  six  temperatures  considered. 

Note  that  according  to  both  the  Schapery  and  Findley  methods  the 
instantaneous  value  of  the  creep  strain,  eQ,  is  a  curve  fitting 
parameter  which  may  be  evaluated  by  our  computer  program  and  is  not  an 
experimental  data  point.  Figure  15  shows  the  variation  of  the 
"theoretical"  transient  creep  strain,  Ae-j ,  with  temperature,  where  Ae-j 
is  defined  as  the  creep  strain  at  30  minutes  U3q)  minus  the  initial 
creep  strain  (eq),  as  determined  by  the  computer  program.  Inspection 
of  Figure  15  will  verify  that  eQ  is  a  curve  fitting  parameter  only, 
inasmuch  as  the  results  for  ae-|  do  not  follow  observed  experimental 
trends  such  as  that  of  Figure  16.  Figure  16  shows  the  variation  of 
the  "approximate  experimental"  transient  creep  strain,  ae,  with 
temperature,  where  As  is  defined  as  the  creep  strain  at  30  minutes 
( ^ 30 )  minus  the  creeP  strain  at  15  seconds  (eq  25)-  Figure  16  indicates 
an  increase  in  transient  creep  strain  with  an  increase  in  temperature 
for  a  given  stress  level.  This  result  is  in  agreement  with  well  known 
behavior  of  viscoelastic  materials. 

It  was  mentioned  in  Chapter  3  that  the  original  computerized 
Schapery  procedure  in  which  the  power  law  exponent,  n,  is  determined 
from  linear  recovery  data  was  found  to  give  inconsistent  results. 
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The  two  stress  invariant  parameters,  n  and  C,  so  determined,  are  plotted 
as  a  function  of  temperature  in  Figures  17  and  18,  respectively.  The 
fluctuation  of  the  parameters  was  so  severe  that  passing  a  curve 
through  these  results  was  not  attempted.  Since  the  goal  of  the  current 
study  was  to  be  able  to  predict  the  Schapery  parameters  for  a  given 
temperature  and  stress  level,  it  is  obvious  that  this  degree  of 
fluctuation  in  the  results  was  intolerable. 

Further  investigation  revealed  even  greater  inconsistencies.  It 
was  discovered  that  the  n-value  increased  as  the  length  of  recovery 
data  increased  for  a  given  length  of  creep.  It  was  also  determined 
that  the  n-value  increased  as  the  length  of  creep  increased  for  a 
given  length  of  recovery.  In  short,  the  power  law  exponent,  n,  is  a 
function  of  both  the  length  of  creep  and  the  length  of  recovery  data 
analyzed.  This  fact  is  quite  distressing,  for  it  implies  that  no 
unique  value  of  the  power  law  exponent,  n,  exists.  This  observation 
could  be  in  fact  consistent  with  the  Schapery  analysis  wherein  the 
reduced  time  parameters  p  and  -p '  might  be  functions  of  time. 

The  author  wishes  to  emphasize  that  the  discussion  of  the  diffi¬ 
culties  with  the  power  law  used  in  the  computerized  version  of  the 
"conventional"  Schapery  procedure  is  relevant  only  for  the  FM-300 
material  system  considered  herein,  and  is  in  no  way  a  statement  of  the 
accuracy  of  the  numerical  algorithm  itself.  The  computerized  "con¬ 
ventional"  Schapery  method  has  been  used  successfully  by  Hi  el  [61] 
for  the  characterization  of  a  934  "neat  epoxy  resin"  used  in  making 
composite  laminates.  He  also  characterized  the  viscoelastic  properties 
of  a  unidirectional  T300/934  graphite/epoxy  laminate  in  the  same  study. 


Stress  invariant  parameter,  C  vs.  temperature 
determined  via  the  'conventional "  Schapery 
procedure,  :  =  452  psi  (3.17  MPa). 
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Modified  Findley  Approach  for  FM-300  Characterization 

In  order  to  stabilize  the  variation  of  the  n-value  with  tempera¬ 
ture,  the  "modified"  Findley  computer  approach  given  in  Chapter  3  was 
employed.  That  is,  the  power  law  exponent  was  determined  from  creep 
rather  than  recovery  data.  The  variation  of  the  parameters,  n,  C,  and 
e  ,  with  temperature,  as  determined  by  this  approach,  are  shown  in 
Figures  19-21,  respectively.  The  variation  of  the  n-value  with  tempera 
ture  of  Figure  19  is  obviously  more  acceptable  than  that  of  Figure  17. 
From  this  point  on,  Figure  19  is  to  be  used  for  determining  the  n-value 
for  a  given  temperature. 

All  of  the  remaining  parameters  necessary  to  evaluate  the  creep 
and  creep  recovery  response  of  FM-300  within  the  temperature  range 
between  room  temperature  and  100°C  (212°F)  are  plotted  in  Figures  22-45 
The  instantaneous  creep  strain,  s  ,  is  plotted  as  a  function  of  stress 
for  each  of  the  six  temperature  levels  considered  in  Figures  22-27. 
Recall  that  £Q  is  a  curve  fitting  parameter  only  and  has  no  experi¬ 
mental  verification.  The  creep  coefficient,  C 1 ,  plotted  as  a  function 

of  stress  for  the  six  temperature  levels,  is  shown  in  Figures  28-33. 
Variation  of  the  vertical  shift  factor,  Ae./g-j,  is  shown  in  Figures 

34-39,  while  the  behavior  of  the  horizontal  shift  factor,  a  ,  is 

a 

detailed  in  Figures  40-45. 

The  nonlinear  parameters,  g  ,  g-j  ,  and  g^ ,  are  shown  plotted  as  a 
function  of  stress  for  various  temperatures  in  Figures  46-57.  These 
figures  are  included  for  inspection  purposes  only,  as  they  are  not 
utilized  in  any  of  the  numerical  calculations  in  the  next  section. 
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Temperature  and  Stress  Situation 


Before  investigating  the  creep  and  creep  recovery  response  of 
FM-300  under  arbitrary  conditions,  it  seems  appropriate  to  pause  for  an 
inspection  of  Figures  58  and  59.  Figure  58  shows  experimental  creep 
strain  data  at  room  temperature  for  several  stress  levels  along  with 
the  associated  fit  of  the  data.  Figure  59  shows  experimental  creep 
recovery  strain  data  at  room  temperature  for  the  same  stress  levels 
along  with  the  Schapery  fit  of  the  data.  4 s  before,  note  that  the 
Schapery  fit  under-predicts  the  recovery  strain  as  time  increases. 

This  brings  us  to  the  point  where  the  most  critical  question  of 
the  current  study  may  be  asked.  Can  the  creep  and  creep  recovery 
response  of  FM-300  be  accurately  predicted  for  any  arbitrary  tempera¬ 
ture  under  the  application  of  any  arbitrary  stress?  Two  arbitrary  creep 
and  creep  recovery  tests  were  performed  in  order  to  answer  this 
question. 

The  first  test  was  performed  at  a  temperature  of  145°F  (62.8°C) 
with  an  applied  stress  of  2260  psi  (15.58  MPa).  The  n-value  cor¬ 
responding  to  this  temperature  was  predicted  from  Figure  19  as  0.126. 

The  instantaneous  strain,  £q ,  was  determined  by  using  linear  interpola¬ 
tion  with  respect  to  temperature  between  Figure  24  and  Figure  25, 
resulting  in  a  value  of  5412.  Similar  linear  interpolations  with 
respect  to  temperature,  between  the  appropriate  figures,  yielded  the 
remaining  necessary  parameters.  The  creep  coefficient,  C' ,  was  found 
to  be  1118,  while  the  values  for  the  horizontal  and  vertical  shifts, 
a  and  Ae,/g,,  were  determined  to  be  equal  to  0.48  and  1414, 
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Experimental  creep  data  (symbols)  and  the  corresponding 
Schapery  fit  of  the  data  (solid  lines)  for  several  stress 
levels  at  T  =  26°C  (79°F). 


Figure  59.  Experimental  creep  recovery  data  (symbols)  and  the 
corresponding  Schapery  fit  of  the  data  (solid  lines) 
for  several  stress  levels  at  T  =  26°C  (79°F). 


respectively.  The  creep  and  creep  recovery  response  was  predicted 
using  equations  (34)  and  (31),  respectively.  Figure  60  shows  a 
comparison  between  the  predicted  creep  strain  via  equation  (34)  and  the 
actual  creep  strain  data,  while  Figure  61  gives  a  comparison  between 
the  predicted  recovery  strain  via  equation  (31)  and  the  actual 
recovery  strain  data.  The  experimental  creep  and  creep  recovery  data 
is  given  in  Table  2. 

The  second  arbitrary  test  was  performed  at  a  temperature  of 
173. 7°F  (78. 7°C)  with  an  applied  stress  of  1808  psi  (12.47  MPa).  The 
five  necessary  parameters  were  evaluated  by  the  same  procedure  as  was 
used  in  the  first  test.  The  n-value,  the  instantaneous  strain,  eQ, 
and  the  creep  coefficient,  C',  were  evaluated  as,  0.165,  4711,  and  776, 
respectively.  The  values  for  the  horizontal  and  vertical  shifts,  a^ 
and  Ae-j/gp  were  evaluated  as,  0.60  and  1166,  respectively.  As  before, 
the  creep  and  creep  recovery  response  was  predicted  using  equations 
(34)  and  (31).  Figure  62  gives  a  comparison  between  the  predicted  creep 
strain  and  the  actual  creep  strain,  while  Figure  63  shows  the  compari¬ 
son  between  the  predicted  recovery  strain  and  the  actual  recovery 
strain.  The  experimental  creep  and  creep  recovery  data  is  given  in 
Table  3. 

Inspection  of  Figures  60-63  reveals  that  the  predicted  creep  and 
creep  recovery  response  was  quite  close  to  the  actual  creep  and  creep 
recovery  response  for  both  of  the  arbitrary  tests  performed.  The  re¬ 
sults  of  these  two  experiments  are  quite  encouraging,  for  they  indicate 
that  the  prediction  of  the  creep  and  creep  recovery  response  of  FM-300 
under  an  arbitrary  temperature  and  stress  situation  is  indeed  possible. 
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Figure  60.  Comparison  between  the  predicted  creep  strain  (solid 
line)  and  the  actual  creep  strain  (symbols)  for 
FM-300,  -  =  2260  psi  (15.58  MPa )  and  T  =  145°F 
(62.8°C) . 


Figure  61.  Comparison  between  the  predicted  creep  recovery 

strain  (solid  line)  and  the  actual  creep  recovery 
strain  (symbols)  for  FM-300,  :  =  2260  psi  (15.58  MPa) 
and  T  =  145°F  (62.8°C). 
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Figure  63.  Comparison  between  the  predicted  creep  recovery  strain 
(solid  line)  and  the  actual  creep  recovery  strain 
(symbols)  for  FM-300,  a  =  1808  psi  (12.47  MPa)  and 
T  =  173. 7°F  (78.7°C) . 


Table  3.  Experimental  Creep  and  Creep  Recovery  Data  for  FM-300  at 
1801  psi  (12.47  MPa)  and  T  =  173. 7°F  (78.7°C). 


Time  (min)  Creep  Strain  ( y/in/in )  Time  (min)  Recovery  Strain 


Material  Parameter  Surfaces 


In  the  preceding  section,  the  desired  material  parameters  were 
determined  by  linear  interpolation  with  respect  to  temperature  between 
the  appropriate  figures.  Determination  of  the  material  parameters  in 
this  manner  implies  the  existence  of  a  material  "parameter  surface"  for 
each  of  the  four  necessary  stress  dependent  material  parameters,  eQ, 

C',  a^,  and  Ae-|/g^.  The  existence  of  the  material  "parameter  surface" 
was  discussed  recently  by  Hi  el  [61],  while  a  similar  concept,  the  lamina 
property  surface  was  discussed  previously  by  Brinson  and  Dillard  [42], 
The  "parameter  surface"  is  a  three  dimensional  illustration  show¬ 
ing  the  variation  of  a  given  material  parameter  as  a  function  of  both 
stress  level  and  temperature.  The  "surface"  for  the  instantaneous 
strain,  e0>  is  shown  in  Figure  64.  Similarly,  surfaces  could  also  be 
drawn  for  the  three  remaining  parameters,  Cf,  a^,  and  Ae-j/g^. 

Possession  of  an  accurate  description  of  these  "parameter  surfaces" 
would  make  the  prediction  of  the  creep  and  creep  recovery  response  of 
FM-300  under  arbitrary  conditions  almost  trivial. 

Prediction  of  Long  Term  Creep  Response  Based  on  Short  Term  Testing 

It  is  well  documented  [61,63,66]  that  the  n-value  derived  via  a 
three  parameter  fit  of  linear  creep  data  increases  as  the  length  of  the 
creep  data  increases,  before  reaching  an  asymptotic  value  at  approxi¬ 
mately  103  minutes.  Although  it  was  known  in  advance  that  predictions 
based  on  the  parameters  which  were  derived  from  only  30  minutes  of 
creep  data  would  under-predict  the  long  term  creep  strain,  it  was  of 
interest  to  verify  this  fact.  Two  long  term  (10“  minutes)  creep  tests 


were  performed  in  order  to  obtain  the  bounds  on  the  accuracy  of  the 
long  term  predictions. 

The  first  long  term  test  was  performed  at  room  temperature  under 
an  applied  stress  of  452  psi  (3.12  MPa).  The  30  minute  creep 
parameters  (eQ  *  610,  C'  *  287,  n  =  0.056)  were  substituted  into  equa¬ 
tion  (34)  and  the  creep  response  through  104  minutes  was  calculated. 

The  comparison  between  the  predicted  creep  response  and  the  actual 
creep  response  is  shown  in  Figure  65.  The  actual  long  term  experi¬ 
mental  creep  data  is  given  in  Table  4.  Note  that  the  actual  creep 
strain  exceeded  the  predicted  creep  strain  by  only  3.0%  at  10,200 
minutes.  - 

The  second  long  term  test  was  performed  at  a  temperature  of  212°F 
(100°C)  under  an  applied  stress  of  1356  psi  (9.35  MPa).  As  before,  the 
30  minute  creep  parameters  (eQ  *  4012,  C'  *  767,  n  =  0.187)  were 
substituted  into  equation  (34)  and  the  creep  response  through  104 
minutes  was  calculated.  The  comparison  between  the  predicted  and  the 
actual  creep  response  is  shown  in  Figure  66.  Table  5  gives  the  actual 
long  term  experimental  creep  data.  In  this  case,  the  actual  creep 
strain  exceeded  the  predicted  creep  strain  by  a  substantial  14.5%  at 
10,246  minutes. 

A  pair  of  interesting  observations  can  be  made  based  on  the  two 
long  term  creep  tests  performed.  First  the  predicted  long  term  creep 
strains  were  always  lower  than  the  actual  creep  strains.  Second,  the 
error  in  predicting  the  long  term  creep  strain  increases  as  the 
temperature  increases. 
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Table  4.  Long  Term  (104  min)  Experimental  Creep  Data  for  FM-300  at 
a  =  452  psi  (3.12  MPa)  and  T  =  79.5°F  (26.4°C). 
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Table  5.  Long  Term  (104  min)  Experimental  Creep  Data  for  FM-300  at 
a  =  1356  psi  (9.35  MPa)  and  T  =  212°F  (100°C). 
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Chapter  6 

SUMMARY  AND  CONCLUSIONS 

The  present  investigation  has  been  concerned  with  the  nonlinear 
viscoelastic  characterization  of  two  commonly  used  structural 
adhesives,  FM-73  and  FM-300,  in  bulk  form.  The  goal  of  the  current 
research  was  to  develop  a  procedure  by  which  the  creep  and  creep 
recovery  response  under  arbitrary  conditions  of  temperature  and  stress 
could  be  predicted  with  reasonable  accuracy. 

A  computerized  Schapery  procedure  was  employed  to  characterize 
FM-73  at  room  temperature.  Several  factors  were  cited  which  may  have 
been  responsible  for  the  inability  to  characterize  this  material  at 
elevated  temperature. 

The  characterization  of  FM-300  via  the  computerized  Schapery 
procedure  utilized  for  FM-73,  was  found  to  yield  inconsistent  results 
for  the  linear  parameters,  n  and  C.  Therefore,  a  "modified"  computer 
procedure  consisting  of  elements  of  both  the  Schapery  and  Findley 
methods  was  utilized.  The  characterization  of  FM-300  via  this  "modi¬ 
fied"  approach  was  successfully  performed  for  six  temperature  levels, 
ranging  from  room  temperature  to  100°C  (212°F). 

It  was  subsequently  shown  that  the  creep  and  creep  recovery 
response  of  FM-300  under  arbitrary  conditions  of  temperature  and  stress 
could  be  predicted  with  reasonable  accuracy  based  on  the  material 
parameters  evaluated  at  each  of  the  six  experimental  temperature  levels. 
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Material  parameters  evaluated  from  30  minutes  of  creep  data  were 
then  used  to  predict  creep  response  over  104  minutes.  It  was  shown 
that  the  predicted  strain  was  always  lower  than  the  actual  strain,  and 
that  this  deviation  between  the  predicted  strain  and  the  actual  strain 
increased  as  the  temperature  increased.  This  under-prediction  of  long 
term  creep  strain  indicates  that  the  value  of  the  power  law  exponent, 
n,  was  too  low.  With  the  method  used  to  determine  n  from  creep  data 
approximately  103  minutes  of  linear  data  are  required  to  determine  its 
value  accurately.  It  might  be  noted  that  Peretz  and  Weitsman  using  a 
similar  method  have  indicated  that  only  fifteen  minutes  of  creep  data 
is  required.  At  this  time  it  is  not  possible  to  reconcile  this 
apparent  contradiction.  However,  regardless  of  the  method  of  evalua¬ 
tion  or  the  duration  of  the  test,  it  is  essential  that  the  n-value  be 
known  accurately  if  long  term  creep  response  is  to  be  predicted 
accurately. 

Although  the  results  of  the  current  research  are  encouraging,  a 
word  of  caution  is  in  order  at  this  point.  The  FM-300  analyzed  in  the 
current  work  was  prepared  according  to  the  cure  cycle  given  in  Chapter 
4.  It  has  been  shown  [67]  that  the  cure  cycle  as  well  as  many  other 
factors  can  have  a  significant  effect  on  the  bulk  tensile  properties 
of  a  structural  adhesive.  Therefore,  results  by  others  on  FM-73  or 
FM-300  processed  differently  and  tested  under  different  conditions 
might  be  at  variance  with  those  of  this  study.  However,  it  is 
appropriate  to  note  that  the  isochronous  stress-strain  curves  for 
t  =  1  min.  obtained  for  FM-73  and  shown  in  Fig.  4  compare  favorably 
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